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1 Introduction 



The title of this symposium, "QCD Corrections and New Physics", provides 
the motivation for this talk: Many potential signatures of new physics require 
a quantitative understanding of Standard Model backgrounds, which often 
can only be adequately provided after QCD corrections to the lowest order 
processes have been computed. The more partons there are in the lowest 
order process, the more important the QCD corrections become. On the other 
hand, even the first corrections, those that are next-to-leading order (NLO) 
in a s , have until recently only been available for processes with relatively few 
external partons, such as single-jet and dijet production at hadron collidersp 
and three-jet production in e + e~ annihilationa 

A major reason for this situation has been the dearth of one-loop QCD 
amplitudes (virtual corrections) for processes with five or more external legs. 
Such amplitudes are in principle straightforward to compute from Feynman 
rules; however, the large number of kinematic variables means that a brute- 
force approach can result in extremely large analytical expressions. Hence 
until recently one- loop amplitudes formed an 'analytical bottleneck' to the 
computation of more Standard Model processes at next-to-leading order. This 
talk will review techniques for systematically disentangling these amplitudes 
into components, primitive amplitudes, which arc sufficiently simple that they 
can be reconstructed essentially entirely from their analytic properties, i.e. 
without resorting to an analytic evaluation of one-loop Feynman diagrams at 
all. 

Some of the techniques discussed here have been used to compute the 
one-loop five-parton amplitudes, ggggg, u qqq'q'g,U and qqgggM which enter 
into NLO studies of the .-three-jet production rate and the substructure of 
jets at hadron colliders&Q More recently, the full set of techniques has led 
to the one-loop helicity amplitudes for a vector bosoa V (which maji.be a 
virtual photon, Z or W) plus four partons: Vqqq'q' B and Vqqgg. The 
helicity-summed versions of the one-loop corrections to qqq' gL^ud ^*qqgg 
were calculated contemporaneously by conventional techniquespJO although 
the analytic formulae were too large to present. After converting-jto a common 
regularization scheme, the two sets of results agree numerically^ 

The Vqqq'q' and Vqqgg amplitudes have hepri used to compute the NLO 
corrections to e + e~ annihilation into four jetsLjli3'E3 The QCD prediction 
for the four-jet rate is almost doubled in going from leading order to NLO, 
bringing it into much better agreement with experiment, for several different 
jet algorithms. Besides the total four-jet rate, jet angular distributions and any 
other infrared-safe event-shape variable whose perturbative expansion begins 
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at order a 2 s can now be computed at NLOBBH The latter results can be 
used to refine experimental measurements of the QCD color factors Cf, Ca 
and Nf . The same amplitudes could also be applied to NLO processes that are 
related by crossing symmetry and coupling constant conversions, such as the 
production of a W or Z plus two jets at hadron colliders, or the deep inelastic 
production of three jets, but this has not yet been carried out. 

Many aspects of the analytical techniques covered here have been pre- 
viously reviewed elsewhereE^SEIlEI] In particular, the recent calculation of 
e + e~ — * V — ► qqgg required a thorough understanding of the spurious as well 
as physical singularities that can occur in six-point amplitudes at one loop; 
such details are beyond the scope of this review. 



2 Primitive Amplitudes and Analytic Properties 

2.1 Traditional Approaches 

In the traditional approach to computing a gauge theory amplitude, one begins 
by drawing all Feynman diagrams. However, individual diagrams are not gauge 
invariant; only certain sums of diagrams are. Thus each diagram contains 
unphysical, gauge- variant pieces, which one would like to avoid calculating 
at all, since such pieces must cancel eventually. Furthermore, the standard 
Feynman rules for the gauge boson self-interactions induce a complicated tangle 
of color and Lorentz indices in each diagram. The three-gluon vertex is 

Vffik,p, q) = f abc (v, P (P - q)» + Vw(q - k) v + v(* - P)p) » (1) 

where f ahc are the SU(3) structure constants, k, p and q the momenta emerging 
from the vertex, and rj^ the Minkowski metric. Each of the three terms in 
cq. (Q) leads to a different routing of the Lorentz indices, and the composition 
of several f abc s similarly leads to different routings of the color indices. 

In a loop amplitude, factors of loop momenta coming from the vertices 
lead to further complications when one attempts to reduce the various loop 
integrals to a basic set of integral functions. The pentagon integral appearing 
in the one-loop five-gluon amplitude includes many terms with five powers of 
loop momenta in the numerator, arising from five vertices of the form ([!]). 
Each such term can amount to a page or more of algebra after reduction of 
the integral. 

The complications of a diagrammatic approach are in marked contrast 
to the simplicity of the final results. For example, in the mid-1980s it was 
found that all QCD tree amplitudes with four or five external partons (and 
many with six or more external partons) could be written as permutation 
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Figure 1: (a) Factorization of a tree amplitude on a particle pole, into the product of 
an (m + l)-point tree amplitude and an (n — m + l)-point tree amplitude, (b) The case 
m = n — 2 (or m = 2) corresponds to two external legs, a and b, becoming collinear. Gray 
blobs represent tree amplitudes, while the hatched blob represents a splitting amplitude. 

sums of simple, single-term expressions (see eq. ([To]))@ Similarly, the one-loop 
four- and fijTCpgluon amplitudes can be represented remarkably compactly (see 
eq. 

2.2 Outline of Analytic Approach 

This contrast suggests that there should be a more direct and manifestly gauge- 
invariant way to calculate loop amplitudes than via Feynman diagrams. Here 
we shall advocate exploiting the two principal analytic properties of loop am- 
plitudes, unitarity and factorization on particle poles, in a 'perturbative_boot- 
strap' approach. In the old S'-matrix bootstrap program of the 1960sp3 one 
attempted to construct S-matrices for the strong interactions using the con- 
straints of analyticity. However, in the absence of a perturbative expansion, 
these nonlinear, coupled constraints proved too difficult to solve exactly, with- 
out imposing artificial constraints of some kind, for example on the number of 
states in the physical spectrum. On the other hand, in the context of pertur- 
bation theory, the constraints of S'-matrix analyticity become recursive in the 
coupling constant g, and therefore tractable. 

For example, an n-point tree amplitude is said to factorize on an interme- 
diate particle pole when the sum of m external momenta satisfies the on-shell 
condition for an intermediate state. The factorization condition for the mass- 
less case is P 2 — > 0, where P = k il + fe^+i + ■ ■ ■ + k i2 . In this limit, depicted 
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in fig. |l](a), the n-point amplitude becomes 



4T(l,2,...,n) 

t° Atree fA A I D\— M _j_ /(tree ( D A • ■ - ■ 

p2 



...,i2, (-P)- X ) ^2 < r - C m+ l(^ A ^2 - 1) , 



where we take all momenta to be outgoing in each amplitude, and there is an 
implicit sum over intermediate polarization states A. (We have suppressed the 
color indices here, as is appropriate for the color-ordered primitive amplitudes 
defined below.) The special case m = n — 2 (or m = 2), shown in fig. |I](b), 
occurs when two of the external momenta become collinear, say k a || k^: 

AT C (- ■ ■ ,a,b, . . .) ^ Spht^V-, b^)AT_\(. . . , P\ . . .) . (3) 

The collinear singularity behaves like the square-root of a pole, which has 
been absorbed into a universal splitting amplitude that replaces the second 
amplitude on the right-hand side of eq. (^). 

The universal structure of these tree-level limits is easy to understand 
heuristically: The spacetime picture of the scattering process reduces to two 
independent scatterings connected by a single intermediate state. It can_also 
be derived from a string-theoretic representation of the tree amplitudesE3 In 
any case, the 'boundary values' (factorization limits) of the n-point amplitude 
contain only lower-point amplitudes or universal functions. Thus the general 
rt-point tree amplitude could be constructed recursively in n, simply by de- 
manding consistency with all factorization limits, provided only that the solu- 
tion to the boundary-value problem is unique. (We will discuss the uniqueness 
issue later.) 

In fact, tree amplitudes can be constructed rather efficiently, and without 
any uniqueness issues, using recursive techniques based on Feynman diagrams Ej'LJ 
However, for loop amplitudes this i s n o t generally the case, and so the factoriza- 
tion properties of loop amplitudesEjo become quite useful. These properties 
are straightforward generalizations of the corresponding tree properties. For 
example, the generic form of the collinear limit of a one-loop amplitude is 

^r 100 ^. . . , a, b, . . .) ^ Split tre A °(« A " , b x ») A\- i° P (- ■ ■ ,P\ ■ ■ ■) (4) 

+ split ^V", 6 Ab X rec i(- • ■ ,P\ ■ ■ ■) , 

where Split 1_loop is a new universal function. 

Even more useful for loop amplitudes are the constraints imposed by uni- 
tarity. After expanding the S'-matrix unitarity relation S^S = 1 in terms of 
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Figure 2: A cut for a one-loop amplitude, with momentum P flowing across the cut. The 
blobs on either side of the cut represent tree amplitudes. 



g, one finds that the cut or absorptive part of a one-loop n-point amplitude is 
given by the two-body phase-space integral of a product of two tree amplitudes, 

4r loo P(l,2,...,n) =i />LIPS(-4,^) A£ e + V(-4r\ii,...,*2,4 2 ) 

P 2 cut J 

x e; + 2((-^)- A2 , 12 + 1, . . . , h - 1, il 1 ) , 

as illustrated in fig. |^. Here the cut channel has momentum squared P 2 , the 
integration is over D-dimensional Lorentz-invariant phase space with with mo- 
menta —i\ and £2 for the intermediate states, and there is again an implicit 
sum over their polarizations Ai,A2- For all the one- loop multi-parton ampli- 
tudes that one might presently contemplate calculating, the tree amplitudes 
required to evaluate each of the cuts are readily available in a compact form. 

2.3 Decomposition into Primitive Amplitudes 

Although the above factorization and unitarity constraints are powerful, they 
are somewhat cumbersome to implement on full multi-parton gauge theory 
amplitudes. It is much simpler to first identify a simpler set of gauge-invariant 
building blocks, called primitive amplitudes™ from which the full amplitudes 
can be built, but which individually have simpler analytic properties. The two 
handles that one can use to pull apart QCD amplitudes are the helicity and 
color quantum numbers of the external quarks and gluons. 

In massless QCD, quark helicity is conserved, and decomposing an ampli- 
tude with respect to it just amounts to inserting a helicity projector ^(1 ±75) 
next to the external spinor. A definite gluon helicity is simplest to implement 
via the spinor helicity formalism£2l which represents gluon polarization vectors 
in terms of Weyl spinors | k }, 

sUk; q ) = ^\n, e - {k . q)= ^ + M k p . ( 6 ) 

mV 1 V2(qk) Ml ; V2[kq] ( ' 
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Here k is the gluon momentum, q is an arbitrary null 'reference momentum' 
which drops out of final gauge-invariant amplitudes, and the spinor inner- 
products are 

(ij) EE (k~\k+) , [ij] EE (k+\kj) , (ij) [J i] = Sij = (h+krf . (7) 

The color structure of SU(N C ) gauge theory amplitudes can be sorted 
into color-ordered components — sums of contributions from Feynman graphs 
where the external states maintain a definite cyclic ordering with respect to 
each other, and where the color-ordered vertices have had the color factors 
removed from themEj The color factors that multiply these components, in the 
color decomposition of the full amplitude, are ordered traces of generators T a 
for the fundamental representation of SU(N C ), Tr(l ...n) = Tr(T ai . . . T a "), 
or strings of the form (T a3 • • • T""),^ 1 if external quarks are also present. For 
example, the color decompositions of the tree-level and one-loop five-gluon 
amplitudes are 

At ee = g 3 Yl TrHl)---cr(5))Ar e (a(l),...,a(5)) (8) 

<tGS 5 /Z 5 



and 



\ 1— loop 
A 5 



N c Tr(a(l) . . . <r(5)) A 5; i((t(1), . . . , a(5)) 



+ T1 M1M2))M<?(3)<?(4:)<J(5)) A 5 ;3(a(l),a(2);a(3),a(4),a(5)) , 

(9) 

where the permutation sums run over all inequivalent traces. j4| roo (l, 2, 3, 4, 5) 
and A 5l i (1, 2, 3, 4, 5) are color-ordered amplitudes, so that when their external 
states are chosen to be of definite helicity, e.g. A 5 ;i(l Al , 2 A2 , 3 A3 , 4 A4 , 5 A5 ), they 
qualify as primitive amplitudes. The same is not quite true of the coefficients 
A 5;3 of the double-trace terms in eq. (^); however, they can be expressed^ as 
sums of permutations of the As ; is. 

The virtue of these decompositions becomes apparent when the explicit 



values of A^ cc are recorded: 

(tree/i ± r>+ Q + 



^ t 5 ree (l ± ,2 + ,3 + ,4+,5 + ) = 0, 
Ar e (l",2-,3+,4+,5+) 



(12> 4 



4 rco (l-,2+,3-,4+,5+) = 



12) (2 3) (3 4) (4 5) (51) ' (10) 
<13) 4 



(12) (2 3) (3 4) (4 5) (51) ' 
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the remaining helicity configurations are related by parity and charge conju- 
gation. 

The zeros in_eq. ( |To| ) are a general consequence of supersymmetric Ward 
identities (SWI)El At tree level, SWI can be applied directly to a non-super- 
symmetric theory like QCD, because no fermions can contribute to a tree-level 
n-gluon amplitude; hence-the 'missing' fermions might as well be considered 
gluinos instead of quarks EHl 

The simple analytic structure of the nonzero terms is due to the color and 
helicity decompositions. Color-ordering implies that factorization singularities 
can only appear in color- adjacent channels. Thus only a subset of the kinematic 
variables — here, {si2, S23, S34, S45, S51} — are 'important' in the sense that 
they contain all poles (and in the loop case, all cuts) in the primitive amplitude. 
The denominator factors in eq. ([To|), (i (i + 1)) = (phase) x ./s^j+i, properly 
capture the square-root behavior of collinear singularities, including a phase 
dependence as ki and fc^+i are rotated around their common collinear axis; both 
of these factors are related to an angular-momentum mismatch of 1 unit in the 
collinear limit. (In other cases, [i (i + 1)] denominator factors can appear, when 
the mismatch has the opposite sign.) 

One-loop primitive amplitudes benefit in precisely the same way as tree 
amplitudes from color and helicity decompositions. In addition, it is possible 
to decompose a one-loop amplitude according to the spins of the particles going 
around the loop, in a 'supersymmetric' way. For example, for a QCD amplitude 
with all external gluons, the internal gluon loop contribution g (and fermion 
loop contribution /) can be rewritten as a supersymmetric contribution plus 
a complex scalar loop s, 

9 = (g + if + Ss) - A(f + s) + s = A N=i - 4 A N=1 + A scalar , 
/=(/ + «)- a = A N=1 - A scalar , ' 

where A N=A represents the contribution of the N = 4 super Yang-Mills mul- 
tiplet, and A N=1 an N — 1 chiral matter supermultipletj.. Related rearrange- 
ments are possible for amplitudes with external quarksl3'E2l) The two super- 
symmetric contributions, A N=A and A N=1 , automatically obey SWI, and have 
other simplifications arising from loop-momentum cancellations. The scalar 
contribution is generally the most algebraically complicated of the three, but 
it is also simpler in some aspects. 

We illustrate this decomposition for A5 ; i(l~, 2 _ , 3 + , 4 + , 5 + ), which is one 
of the two one-loop primitive amplitudes required to calculate the 'gluonic' 
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NLO corrections to pp — > 3 jets. Its components according to (O) are 



.AT=4 



i=i 



. ^_ /ttrcc 
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1 








['»( 




"2 





' 1,2 



S23 



hi 



hi 



Sj + l,j+2 
2 



hi 



-Sj+2,j-2 
Sj-2,j-l 



7T" 



-ssi 



+ 



icr (1 2) 2 ((2 3) [3 4] (4 1) + (2 4) [4 5] (5 1» ln 



(2 3) (3 4) (4 5) (51) 



S51 - s 2 3 



A' 



calar 



-A N=1 + — 



2 Al iee 



[3 4] (4 1) (2 4) [4 5] ((2 3) [3 4] (4 1) + (2 4) [4 5] (5 1) ) 



In I 



( -S23 ] —if £23 _ s 51 ] 

\_ — SB1 / 2 y S51 S23 / 



(3 4) (4 5) 



(3 5) [3 5] J 



(12) [3 5f 



[12] [2 3] (3 4) (4 5) [51] [2 3] (3 4) (4 5) [5 1] 



(s 5 i - S23) 3 
(12) [3 4] (41) (2 4) [4 5] 



S23 (3 4) (4 5) S51 

(12) 

where c r = T(l + e)r 2 (l - e)/((47r) 2 ~T(l - 2e)). Since the three components 
have quite different analytic structure, the rearrangement ( pi] ) is a natural one. 
As expected, the A N=A is the simplest component, followed by A N=1 . While 
^scalar j g most complicated piece, it has (for example) no S34 cut, while 
the full gluon loop does. 

Although we do not have space to demonstrate it explicitly here, evaluation 
of the unitarity cuts, eq. (|^), for an amplitude such as eq. dl2| ) is quite simple. 
The principle reasons are: 

(1) The tree amplitudes are fully 'processed', in terms of gauge cancellations, 
etc., before they are fed into the cut. 

(2) SWI generally apply to the tree amplitudes, simplifying their form, even if 
they do not directly apply to the one-loop amplitudes. 

(3) On-shell conditions for the intermediate legs, i\ = l\ — 0, can be used 
repeatedly to simplify the evaluation. 

The only catch to a unitarity-based approach is that loop amplitudes can 
have rational-function terms, containing no logarithms or dilogarithms, which 
are not detectable in any cut (when the cut is evaluated in four-dimensions). 
In eq. (p^|), such terms appear only in the non-supersymmetric component 
^scalar _ ^phe rational- function terms in the supersymmetxic components can 
be shown quite generally to be linked to the logarithmsEj) Fortunately these 
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terms can be determined using the factorization properties of loop amplitudes, 
such as eq. (^). One can construct an ansatz for the rational- function terms, 
and then proceed channel by channel to add terms to the ansatz so that it 
correctly reproduces the desired singular behavior in each channel, using fac- 
torization information provided by lower-point amplitudes. This procedure 
could fail if the factorization properties did not uniquely fix the amplitude. 
Although we have no formal uniqueness proof, empirically the procedure al- 
ways works for six-point amplitudes. We can verify that it has worked by 
numerical evaluating conventional Feynman diagrams at one or more random 
kinematic points. (At the five-point level, there is a special term that provides 
a counter-example to uniqueness!^! although in practice its coefficient could 
also be fixed numerically.) 

The above strategy of calculating primitive amplitudes via their unitarity 
cuts, with the rational-function terms fixed via factorization, has been success- 
fully applied to the one-loop e + e~ — * 7, Z — + 4 parton amplitudes. Most of 
the previously calculated one-loop four- and five-parton amplitudes reappear 
in this calculation as the boundary-value information for the rational-function 
terms, thus emphasizing the bootstrap nature of the approach. 

3 Phenomenological Applications 

3.1 pp — > 3 jets 

To date there have been two phenomenological applications of the one-loop 
helicity amplitudes whose calculation was outlined above. First, two groups 
have undertaken to compute NLO corrections to three-jet production in pp 
collisionsQiJ Production of three jets is one of the dominant processes at large 
transverse momentum at the Tevatron, after thex>ne-jet inclusive and two-jet 
rates which have already been computed at NLOEI By measuring the three-jet 
to two-jet ratio, once the former has been computed at NLO, one should be able 
to extract a relatively precise value for a s {Q 2 ), at the highest experimentally 
accessible values of Q 2 . 

Although all the matrix elements required for the full NLO pp — + 3 jets 
calculation are available, the extreme computational demands of the numeri- 
cal integrals encountered in the calculation have led both groups to perform 
a 'gluonic approximation' to the true result, in which the quark distribution 
functions for the proton are dropped, as_weU as final-state quarks. Thus only 
the o nfirjQQD five-gluon matrix elemental and tree-level six-gluon matrix ele- 
mentsooo have to be evaluated. Although this 'approximation' is not ex- 
pected to be particularly good quantitatively, it provides some useful qualita- 
tive lessons. 
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Figure 3: Dependence of the gluonic pp — > 3 jet cross-section on the jet resolution parameter 
cfcut for the fcy clustering algorithmic The gray and black bands result from varying the 
renormalization/factorization scale ji over the range d cu t/2 < fi < 2d cu t for the Born and 
NLO results, respectively. 

TrocsanyiE] used an iterative fcx clustering algorithmB adapted from e + e~ 
annihilation, with a variable jet resolution parameter d cu t which controls the 
overall hardness of the three-jet event. The d cu t-dependence of the leading- 
order (Born) and NLO results are shown in fig. ||. For this jet algorithm, the 
Born and NLO cross-section have very similar ei cu t-dependence. However, the 
NLO result has (as expected) a much reduced dependence on the (unphysical) 
renormalization scale fi, which bodes well for the accuracy of the full NLO 
three-jet rate once quark contributions are included. 

Kilgore and Gielei studied pp — > 3 jets for the kx algorithm as well as for 
three algorithms of the cone type more commonly used at hadron colliders: 

1) a 'fixed cone' algorithm, used by UA2J23 

2) an 'iterative cone' algorithm, used by CDF and D0,EZl 

3) the 'EKS' algorithm, used in NLO one- and two-jet inclusive calculationsEa 
Their NLO results for all but the iterative cone algorithm are shown in fig. [|. 
The dependence of the three-jet cross-section on the transverse energy of the 
leading jet, changes considerably between leading and next-to-leading 
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Figure 4: The NLO gluonic pp — > 3 jet cross-section as a function of the transverse energy 
of the leading jet, for three different jet algorithms and fi = 100 GeV, from ref. ^. See 

text for descriptions of the algorithms. The Born (LO) result is also shown for comparison. 

order; it also varies with the jet algorithm. 

There is no plot for the iterative cone algorithm — the one presently used 
by CDF and DO — because Kilgore and Giele found it to be infrared unsafe for 
the three-jet cross-section. (See also ref. |39|.) They traced this problem to final 
states with three hard partons, two of which are separated by slightly more 
than the cone size R. If a soft gluon is added to this configuration somewhere 
between these two hard partons, its only effect in the other jet algorithms will 
be to shift one of the three jet axes slightly. But in the iterative cone algorithm 
this shift will (in a subsequent step) trigger the merging of two of the three jets 
into one, thus altering the three-jet rate. (In fact, CDF and DO had already 
introduced an additional jet separation cut when comparing their multi-jet 
data to leading order predictions, which effectively removed this problem, at 
least at NLO.) The full pp — > 3 jet NLO cross-section still awaits the inclusion 
of quarks; however, the infrared dangers of the iterative cone algorithm already 
constitute a useful lesson from the gluon-only study. 
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3.2 e+e — > 4 jeis 

A second phenomenological application has been in e + e~ annihilation. Here 
four-jet production, or more generally infrared-safe event-shape variables whose 
perturbative expansion begins at order oq, can now be studied to NLO accu- 
racy. Such observables are sensitive to the non-Abelian self-coupling of gluons 
and to the production of hypothetical colored, electrically neutral particles 
such as light gluinos. Also, at LEP2 energies e + e~ — > 7, Z — » 4 jets forms a 
significant background to W pair production when both Ws decay hadroni- 
cally. 

To date, two independent numerical programs for NLO corrections to 
generic order-af event shapes have been constructed, MENLO_PARCe2i and 
DEBRECENPj implementing the one-loop matrix elements for e + e~ — * 7, Z-.f^, 
4 partons, and the tree-level matrix elements for e + e~ — * 7, Z — + 5 partonsrHr 2 ! 
Althoiiffik-Lhe two programs use different generalizations of the subtraction 
methodocj for carrying out integrals over the singular five-parton phase-space, 

itistical errors for the Monte Carlo integration. 
3 for the four-jet fraction i?4 using the Durham 



their results agree to within s 
In fig. ||, QCD predictions 
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jet algorithnE§E3 are compared with data at the Z pole from ALEPHeZI and 
SLDpj as a function of the jet resolution parameter ?/ cut . It has long been 
known that the Born-level QCD prediction for the Durham (and most other 
e + e _ jet) algorithms, when evaluated at a renormalization scale fi = \fs, 
under-predicts the measured four-jet fraction by roughly a factor of two. The 
figure shows that addition of the NLO correction brings the prediction into 
much better agreement with the data, within roughly 10%. 

For small values of y C ut, where the four-jet fraction is large, a kinematic 
singularity develops in the perturbative expansion, and the true expansion pa- 
rameter becomes a s L 2 , where L — ln(l/y cut ). The Durham algorithm has the 
virtue of being resummable — the leading and next-to-leading logarithms at 
each order in a s (terms of the form a™L 2 ™ and a™L 2n_1 ) can be calculated. In 
fig. H we also show the result of carrying out this resummation, and 'matching' 
the result to the fixed-order (NLO) calculation, which amounts to subtracting 
out common terms in the two expansions. This one-loop plus resummed result 
does improve the agreement with data still further (for fi = Mz) at small y cu t- 

Although it has been suppressed in fig. ^|, the NLO prediction for R± 
still has a reasonably large uncertainty stemming from the truncation of the 
perturbation series after just two terms, and the large size of the second term. 
(This uncertainty is perhaps 20-30% if one estimates from the renormalization- 
scale dependence of the NLO result.) This feature limits the utility of R4 in 
precision tests of QCD. On the other hand, normalized distributions for various 
angles between the four jets, tend to have quite small, and hence reliable, NLO 
corrections, of order 1-3%L3'E3 



4 NNLO Predictions for Jet Physics? 

Uncertainties from higher-order terms still plague NLO calculations, though 
not as severely as at Born level. For example, a s (Af§) has been measured 
in e + e~ annihilation at the Z pole by using a dozen or more event-shape 
observables Oi whose perturbative expansion begins at order a s (such as the 
three-jet fraction R3): 




The extracted values of a s (Af|), obtained by fitting to the 0(a 2 s ) truncation 
of eq. (filf), depend significantly on the value assumed for fi (an unphysical 
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parameter), and on the observable used, leading to a theory uncertainty that 
dominates the overall errorpa 

a s (M|) = 0.121 ±0.002(exp.)± 0.005(theory). (14) 

Calculation of the next-to-next-to-leading order (NNLO) 0{a s s ) terms in eq. ( p"g| ) 
would clearly improve this situation. 

To date, NNLO predictions in QCD are available for only a limited num- 
ber of rather, inclusive quantities, such as the total cross-section in e + e~ 
annihilations (Because such observables have a perturbative expansion of the 
form 1 + a s /Tr + ■ ■ •, a competitive a s extraction requires a very high precision 
measurement of the observable.) NNLO results for jet physics are still a ways 
off; however, some of the ingredients are starting to be attacked. For the above 
observables Oi, three major classes of matrix elements are required: 

(1) tree-level e + e _ — * 7, Z — ► 5 partons, 

(2) one-loop e + e~ — > 7, Z — > 4 partons, 

(3) two-loop e + e~ — > 7, Z — > 3 partons. 

As we have seen, the first two classes are now available. There will also be a 
considerable amount of work required in order to reliably integrate the second, 
and particularly the first, class of contributions over singular ^prners of phase 
space. Some steps have recently been taken along these linesEil 

As for the two-loop matrix elements, an encouraging sign is that a complete 
two-loop four-gluon scattering amplitude has been computed for the first time, 
using generalizations of the analytic tools described above, and the result is 
quite compact^ There are two caveats: the theory is not QCD, but N = 4 
supersymmetric SU (N c ) Yang-Mills theory; and the answer has been expressed 
in terms of scalar integrals, but those integrals have not yet been performed in 
closed form. The result for the leading-color part of the amplitude is simply 

A N=4, 2-loop, L.C. = _ N 2 stA tree q + f ^ ^ ^ 

where I±(s, t) is the two- loop planar double-box scalar integral, 
•d D pd D q 1 



X 4 (s,t) 



{2ttY d P 2(p - ki) 2 (p - k t - fc 2 ) 2 (p + <z) V(« " h) 2 (q ~ h - fc 4 ) 2 

(16) 

Whether two-loop results can be obtained in the same way for QCD remains 
to be seen. 

5 Conclusions 

In charting the progress of particle theory from the 5-matrix days of the 1960s, 
to the triumph of the Standard Model in the 1970s, to the recognition of its 
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probable role as an effective field theory in some grander theory, WeinbergLJ 
remarked: 

The justification of any particular effective field theory is that it is simply 
the most general possible theory that satisfies the axioms of analyticity, 
unitarity, and cluster decomposition along with the relevant symmetry 
principles, so in a way our use today of effective field theories is the 
ultimate revenge of S-matrix theory . . . Quantum field theory is nothing 
but S'-matrix theory made practical. 
In some sense, the techniques reviewed here represent a further revenge of 
S'-matrix theory at the computational level. Although the purely S-matrix- 
based bootstrap program of the 1960s proved intractable, the lesson here is 
that a perturbative bootstrap program can succeed. Once helicity, color and 
supersymmetry decompositions are used to break up loop amplitudes into com- 
ponents with simpler analytic structure, this approach becomes a particularly 
efficient way to compute. As a practical result, improved QCD predictions for 
various multi-jet processes are now emerging. 
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